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INTRODUCTION 

Previous studies of the radiation of sound from vibrating cylindrical 
beams have been concerned with radiation from resonant modes of these finite 
beams, with specified end boundary conditions. These include studies by 
Yousri and Fahy (ref. 1) and Kuhn and Morfey (ref. 2). The radiation efficien- 
cy or radiation loss factor determined in this manner represents the con- 
tribution of a single mode to the radiation at a given frequency. More 
recently, Yousri and Fahy have presented (ref. 3) a more general derivation 
which shows that the radiation efficiency for a cylindrical beam is a sum- 
mation of terms that represent contributions from various modes. This summing 
over modes is necessary whenever more than one mode is excited in the fre- 
quency band of Interest. 

The results of the present study are given in a form which depends only 
on the frequency of the beam vibrations and the physical characteristics of 
the beam and its surroundings. A statistical consideration of random beam 
vibrations allows this result to be independent of the boundary conditions at 
the ends of. the beam. The acoustic power radiated by the beam can be deter- 
mined from a knowledge of the frequency band vibration data without a 
knowledge of the individual modal vibration amplitudes. 

A practical example of the usefulness of this technique is provided by 
the application of the theoretical calculations to the prediction of the 
octave band acoustic power output of the picking sticks of an automatic 
textile loom. Calculations are made of the expected octave band sound pressure 
levels based on measured acceleration data. These theoretical levels are 
subsequently compared with actual sound pressure level measurements of loom 
noise . 


THEORY 

A beam of finite length is modelled as a cylinder of infinite length 
situated on the z axis (see Fig. 1). The transverse velocity of the beam 
vibrations is assumed to be zero except on the segment of the cylinder which 
lies between the points z-l/2 and z-- % H . On this vibrating segment 
the transverse velocity consists of x and y components, v and v , which 
are expanded in a Fourier series of arbitrary fundamental time^period 


1071 



( 1 ) 


ns -ao 

The individual components of the Fourier series are represented as the super- 
position of two traveling waves moving in opposite directions on the beam 
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where k^ is the wave number of the beam vibrations, and represents the 
relative phase between the two traveling waves. Furthemore, v (z) can be 
written as 


v^(z) 


a. 


) e dx 


and the radial or normal velocity is 


Vj^(0,z) = (cosQ t + sine 


v( a ) e^*^ ^ do 
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The acoustic pressure due to these transverse vibrations is found by 
applying the acoustic boundary condition at the surface of the cylinder to the 
solution of the linear acoustic wave equation in cylindrical coordinates. The 
partial differential equation and the appropriate boundary condition are 


^ 4. ^ is 4. J: 4. 
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The solution to Eq. 5 can be written as a general linear combination of the 
seperable solutions (see ref. 4). The solution for outgoing waves of fixed 
frequency ii> ^ n OJe can be written 
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P(r,e,z) \ [Aj.a )cos(m0) + B^<a )sin(me) do( 
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where k =(k^- a^) , k=co /c and is the Hankel function of the first 

kind and^of order m . Applying the™boundary condition of Eq. 6 to the above 


expression, one finds that 



P(r,0,z) = (cos0 t + sin0 J) • 

1 Z( a,r)v(a )e da 
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The time average acoustic power radiated by the beam at a given frequency 
is found by taking the time average of the integral over the surface of the 
cylinder of the product of the acoustic pressure at the surface of the cylinder 
P(r^,e,z) and the normal velocity v (0,z) 


W 



P(r^,e,z) v*(0,z) r^d0 dz 


( 10 ) 


where Re denotes the real part, and * indicates the complex conjugate. The 
ensemble average acoustic power radiated in a given frequency band is found by 
taking the ensemble average of a sum over the frequencies within the band of 
the result of Eq. 10. The ensemble average is performed assuming all relative 
phases to be equally probable. 


The result of the 0 integration in Eq. 10 is TT , and the z integra- 
tion produces a dirac delta function 2TT6(a-a‘') , This allows one of the a 
integrations to be performed by inspection, yielding 

W = If^ro \ Re {Z^(a)} (|v^(a)|^ + |v (a)|^) da (11) 
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This expression for the acoustic power radiated by a finite cylinder can 
be simplified by applying some of the properties of Hankel functions. For 
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k < a the |rgiimgnt for the Hankel functions is imaginary and Re{Z (a)} is 
zero. For k > a the argument is real and the Wronskian relation ?or the 
Hankel functions gives (ref. 5) 
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Referring back to Eqs. 2 and 3, one can write v(a) as the inverse 
transform of v(z) , which when evaluated gives 
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The cos'll term goes to zero in taking the average assuming all ip to be 
equally probable. It can be shown, by taking the time average of v«v , that 
the mean square velocity in a given frequency band is equal to a Sum over the 
frequencies in the band of the magnitude squared of v^ 
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assuming x and y vibrations to be uncorrelated. 

All these results are incorporated into Eq. 11 with the result that the 
acoustic power radiated in a given frequency band is 
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where all frequency terms are evaluated at the center frequency of the band. 

This result can be written in the form of a radiation efficiency for the 
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beam. For a cylinder the radiation efficiency is 


a = 
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With the change of variables 3 = oi/k and the above expression, one can write 
the radiation efficiency for a cylindrical beam 
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where e » k^/k and a symmetry of the integral was used to simplify the 
expression. 

RESULTS 

In the low frequency or small radius limit , kr is small . An approx- 
imation may be used in place of the Hankel function in order to simplify the 
integral in Eq. 19. The first term in the series expansion for the derivative 
of the Hankel function gives 
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With this approximation and the change of variables u = (e+3)kJl/2 , one gets 


a = 


(kr^) 

o 


3 f 


(e-l) 


ki 


r, /2u ^^2, sin u 

[1 - (jj - e) 1 — 2- du 

U 


( 21 ) 


Well below the coincidence frequency (k, =k) » u in the denominator can 
be approximated as u2 = (kJL/2)^ , The result of this approximation is 


a = I (kr„)^ (^) (f-)^ [1 - 
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This expression is directly comparable to the results of Kuhn and Morfey (ref. 
2) in their low frequency approximation of the Yousri and Fahy (ref. 1) ex- 
pression for the radiation efficiency of a simply supported beam. The radia- 
tion efficiency for the simply supported case gives a result which is exactly 
twice that given by Eq. 2 2 when evaluated at the same resonance frequency. The 
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reason for this difference lies in the fact that the simply supported boundary 
condition Implies a specific phase ^ in Eq. 2. This phase happens to be one 
which maximizes the radiation efficiency. Other possible phase relationships 
result in lower values for the radiation efficiency. It will be shown later 
that these differences disappear at higher frequencies. 

At High frequencies the major contribution to the integral in Eq. 19 
comes from the vicinity of 3 “ -e • Expanding the Hankel function term in a 
Taylor series about that point and keeping only the first term one gets 

„ . _2 I (23) 

■ "‘"o (l-e2)|g(tt'(^2’k, )|2 

'1 O' 


This is the radiation efficiency for a cylinder of infinite length. If 
(l-e2)^'^kr » 1 , then the derivative of the Hankel function may be expressed 
in terms of its asymptotic limit and 


0 = 
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(l-e^)(kr^)^+l 


(24) 


For purposes of a numerical evaluation of Eq. 19, the radiation efficiency 
is considered as a function of three independent variables only one of which is 
frequency dependent. These are kr^ , and 


e(kr^) 


1/2 ^ ,a;l/4 
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where m is the mass per unit length of the beam and B is the bending 
stiffness. This last parameter is chosen so as to eliminate .the frequency 
dependence of e = k, /k . Figs. 2 and 3 show the results of a numerical 
calculation of the radiation efficiency. Each graph shows 0 as a function 
of kr for several values of e(kr )^'2 for one value of ^/r^. A 

valid comparison with the modal approach can be made by leaving the phase angle 
in Eq. 15 and continuing the derivation of the radiation efficiency. The result 
is 
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Fig. 4 gives this comparison for the two extreme cases, where the cosip term 
is either always positive or always negative at the modal resonance frequencies. 
It is clear from the figure that the contribution from the second term to the 


1076 



integrand in Eq, 26 diminishes at higher frequencies, as the two extreme cases 
converge towards the average. 

APPLICATION TO TEXTILE LOOM PICKING STICKS 

In order to apply these results to the picking sticks of an automatic 
textile loom, it is necessary to generalize them to allow for different 
radiation efficiencies and perhaps different effective radii for the x and y 
vibrations. This is necessary to account for the fact that the picking stick 
is more nearly rectangular in cross section than circular. With this in mind 
one can write the acoustic power as 
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Values for the parameters were chosen to correspond to the characteristics 
of the picking sticks. The numerical technique used to evaluate O and 
measured vibration data were used to determine the octave band acoustic power 
output for each of the two picking sticks on a loom. Fig. 5 shows a graph of 
the resulting power levels. 

From these power levels, octave band sound pressure levels were calcula- 
ted at a reference point one meter from the front of the loom, assuming 
symmetric cylindrical spreading, 

SPL = X0 1og^„^ <p2> - ^ ["left “right' 

ref 

—5 2 

where R is the distance to the reference point and 2 X 10 N/m is 

the reference pressure. A comparison of this predicted sound pressure level 
with sound pressure levels actually measured at this position is shown in 
Fig. 6. This graph shows good agreement between theoretical and experimental 
results for frequencies above 125 Hz. The agreement is particularly close 
in the range of frequencies which have the highest Sound pressure levels. Both 
curves in the figure represent an overall A-weighted level of 94 dBA. It is 
clear however that the low frequency predicted result falls short of the 
measured values. It is thought that there may be other noise sources on the 
loom which contribute to the higher levels at these low frequencies. 
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Figure 5.- Theoretical acoustic power output. 
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Figure 6.- Theoretical and experimental sound pressure levels. 
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